Quantum Phase Diagram and Excitations for the One-Dimensional S — 1 Heisenberg 

Antiferromagnet with Single-Ion Anisotropy 
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We investigate the zero-temperature phase diagram of the one- dimensional 5=1 Heisenberg 
antiferromagnet with single-ion anisotropy. By employing high-order series expansions and quantum 
Monte Carlo simulations we obtain accurate estimates for the critical points separating different 
phases in the quantum phase diagram. Additionally, excitation spectra and gaps are obtained. 
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I. INTRODUCTION 

Interest in one-dimensional 5 = 1 antiferromagnets is 
long-standing and can be traced back to the original work 
by Haldane* 1 ^ By analyzing the presence of topologi- 
cal terms in effective field-theories for one-dimensional 
antiferromagnets he conjectured that integer-spin chains 
display a ground-state with exponentially decaying spin- 
spin correlations and a gapped excitation spectrum, 
properties markedly different from those displayed by the 
exactly solvable S = 1/2 chain. Despite early contro- 
versy, this so called "Haldane conjecture" is now sup- 
ported by solid numerical 3 - and experimental 4 ^ evidence 
(for a review of earlier results, see Ref. @). 

We are interested in the anisotropic S = 1 antiferro- 
magnctic chain described by the Hamiltonian 



(1) 



The single-ion anisotropy term proportional to D is rel- 
evant in accounting for the magnetic properties of a 
number of compounds: CsNiCl3 (weak axial anisotropy, 
Rcfs. and !), NENP [NifC 2 H 8 N 2 ) 2 N0 2 (C10 4 )] 
(weak axial anisotropy, Ref. (8|), CsFeBr 3 , NENC 
[Ni(C 2 H 8 N 2 ) 2 Ni(CN) 4 ] and DTN [NiCl 2 - 4SC(NH 2 ) 2 ] 
(strong planar anisotropy, Refs. H, Qj] and [ll| re- 
spectively); DTN is particularly interesting due to 
Bose-Einstein condensation of spin excitations under 
a magnetic fieldi 11 ' 12 The model defined by Eq. ((T|) 
is also very appealing from a purely theoretical per- 
spective: the single- ion anisotropy D/J controls the 
magnitude of quantum fluctuations and stabilizes dif- 
ferent phases that, along with the quantum crit- 
ical points separating them, have been extensively 
investigated . 13 - 14,15 - 16,17 - 18,19 - 20,21 - 22,23 - 24,25 

For large positive values of the single-ion anisotropy 
D/J the system is in an Ising-like Neel phase character- 
ized by finite staggered magnetization at zero tempera- 
ture T = and displaying a gap for holon-like excitations, 
that we will discuss later in Sec. IIII B[ Comparatively 
less well understood is the so-called large-D phase for 
large negative values of D/J, first investigated in detail 
by Papanicolaou and Spathis^ The ground-state in the 
large-D phase smoothly evolves from the ground-state 



for D/J — > — oo, simply given by the tensor-product of 
Sf = states at each site. The lowest-energy excita- 
tions in this phase reside in the — ±1 sector and were 
termed excitons and anti-excitons in Ref. [l6l . where the 
existence of bound-states was also verified. 

Physical intuition was gained about the intermediate 
phase observed for small values of \D\/J, including the 
isotropic point D/J = analyzed by Haldane^ by 
studying an extended S = 1 model with biquadratic in- 
teractions 



Waklt — J ^ \Si ■ Si + \ + f3(Si ■ §i+iY 



(2) 



Affleck, Kennedy, Lieb and Tasaki (AKLT)2& showed 
that this Hamiltonian is exactly solvable at (3 = 1/3, 
where it displays a simple valence-bond solid (VBS) 
ground-state with gapped excitations, den Nijs and 
Rommelse^ 7 . later suggested that this VBS state can be 
interpreted as a "fluid" with positional disorder and long- 
range antifcrromagnctic order, characterized by a finite 
expectation value for the string-operator, 



0.(r) 



■- Sq exp 




(3) 



in the limit r — > oo, associated with the breaking of a hid- 
den Z 2 x Z 2 symmetry^ Since the ground-state at = 
has been shown 3 - to exhibit long-ranged string correla- 
tions and is adiabatically connected to the ground-state 
at j3 = 1/3 (sec e.g. Ref. |29( ). one concludes that the 
Haldane phase has a VBS character. We remark that in- 
terest in Haldane-like phases exhibiting long-range string 
correlations has been renewed and proposals for investi- 
gating string-order in cold atomic systems have recently 
been made i 30 i 31 

In this paper, we are primarily interested in improv- 
ing on previous estimates for the location of the largc- 
D — Haldane (DH) and the Haldane— Neel (HN) critical 
points. Improved results may be used as benchmarks for 
further tests on the applicability of quantum information 
tools to detect quantum phase transitions in the model 
Eq. |T]) . 24 ' 25 Additionally, we obtain results for the exci- 
tation spectra and gaps in the large-D and Neel phases 
that may be of experimental relevance. We also show 



2 



that a series expansion method previously applied for 
calculating the Haldane gap22 yields an incorrect result, 
and we propose a different approach for future work. 



II. METHODS 

We have investigated the 5 = 1 antiferromagnetic 
chain described by the Hamiltonian Eq. |T]) by combining 
high-order series expansions and quantum Monte Carlo 
(QMC) simulations. Technical details are described in 
this Section. 



A. Series Expansions 
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Numerical linkcd-clustcr expansions have been extcn 
sively employed in the investigation of quantum magnets 
(for detailed accounts the reader is referred to Refs. 
and 
and 
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341 ; a closely related method is discussed in Refs. 
361 ) and have been recently derived by some of u: 
for the square lattice version of the Hamiltonian Eq. |T]). 
Among the method's many appealing features we high- 
light its applicability to the study of excitations and dy- 
namical responses, a notoriously difficult task within al- 
ternative numerical approaches, following the procedure 
originally devised by Gelfandi^ 

Basically, the linked-cluster method relics on a suitable 
decomposition of the lattice Hamiltonian under investi- 
gation, 



H = H + XV 



(4) 



It is assumed that the ground-state of the unperturbed 
Hamiltonian Tlo is known and can be written as a ten- 
sor product of local states. One proceeds by deriving 
a standard Rayleigh-Schrodinger perturbative expansion 
for connected clusters comprised of an increasingly large 
number of sites and, at each step, one subtracts con- 
tributions from embedded subclusters containing fewer 
sites. In this way, long series for ground-state quantities 
and excited states can be exactly calculated up to a cer- 
tain order in the expansion parameter A, that are then 
analyzed by means of any suitable extrapolation tech- 
nique (in this paper we adopt a standard Pade analysis). 
In what follows, we present the different expansions em- 
ployed in the present study. 



FIG. 1: (Color online) Ground-state energy density for the 
model described by Eq. ([TJ as a function of the single-ion 
anisotropy D/J. QMC results are shown as filled circles 
and have been obtained for a lattice with L = 48 sites. 
Empty symbols correspond to results from the series expan- 
sions: large-D expansion (LDE, diamonds), Neel expansion 
(NE, squares) and staggered-field expansion (SFE, triangles). 
The vertical dashed lines highlight our best estimates for 
the location of the critical points: Dq U /J — —0.971(5) for 
the point separating the large-D and Haldane phases, and 
D§ N /J = 0.316(2) for the HN phase transition (see main 
text). Unless visible, error bars are smaller than depicted 
symbols. 



the superexchange terms in Eq. ([I]) being treated as the 
perturbation y LDE 



V 



LDE 



QZ QZ 



i+1 



(6) 



The (non-degenerate) ground-state of 7ig DE is simply 
given by tensor products of Sf = states for all sites 
i on the chain and, in what follows, we refer to this ex- 
pansion as the large-D expansion (LDE) . Results for the 
ground-state energy density obtained from an LDE series 
comprising terms of up to A 30 are shown in Fig. [T] (open 
diamonds). Shorter series have been obtained for the ex- 
cited states (up to A 14 ), the results being presented in 

Sec. Ens 



2. Neel Expansion 



1. Large-D Expansion 

For investigating the phase stabilized for large nega- 
tive values of D / J we consider the single- ion term as our 
unperturbed Hamiltonian 



A suitable linked-cluster expansion for large positive 
values of D/J is obtained by choosing the unperturbed 
Hamiltonian 

nr = j^s^s^-d^sD 2 , (7) 



(5) and starting from a perfectly ordered Neel unperturbed 
ground-state. The remaining fluctuation terms in Eq. ((T|) 
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are our perturbation V 



NE 



(8) 



and we refer to the resulting expansion as the iVeeZ ex- 
pansion (NE). The ground-state energy and the stag- 
gered magnetization in the Neel phase have been cal- 
culated from a NE comprising terms of to A 18 (results 
for the ground-state energy are depicted as open squares 
in Fig. [l}. Scries for the single-particle excitations have 
been obtained from expansions up to A ; the results are 
discussed in Sec. IIIIBl 



3. Staggered- Field Expansion 
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We have also considered a series expansion originally 
devised by Singh,— here dubbed as the staggered- field 
expansion (SFE), obtained by considering a Hamiltonian 
including an artificial staggered-field term 



n 



SFE 



QZ QZ 



2 {^t^i+l + &i ^i+l) 



(9) 



(Terms proportional to A are our perturbation V and 
the remaining ones define the unperturbed Hamiltonian 
Hq FE .) The unperturbed ground-state is a perfectly or- 
dered Neel state, as in the case of the NE, and the Hamil- 
tonian Eq. (TT]) is recovered in the limit A = 1, where the 
last term in Eq. ((9]) vanishes. This staggered-field term is 
included in order to prevent the condensation of solitons 
before A = 1, as discussed in Ref. l32|, allowing one to 
obtain a seemingly precise estimate for the Haldane gap 
at D/J = 0. Results for the ground-state energy per site 
obtained from an SFE including terms up to A 18 are rep- 
resented by open triangles in Fig. [T] it is clear from this 
plot that the SFE systematically underestimates Eq/L 
towards the DH critical point and we will address this 
issue later in Sec. IIII C\ along with our results for the 
gap at D/J — 0. 



B. Quantum Monte Carlo 



FIG. 2: (Color online) (a) Spin-stiffness ps as a function of the 
single-ion anisotropy close to the critical point separating the 
large- D and Haldane phases (lines are only guides to the eye). 
Results have been obtained from QMC simulations for values 
of L ranging from 24 to 48 and error bars are smaller than 
depicted symbols, (b) Data collapse is achieved for Dq 11 /J = 
-0.971(5) and v = 1.4(1) (see main text). 



are converged, within statistical errors, to their ground- 
state expectation values for this value of T, as verified by 
running preliminary QMC simulations for various tem- 
peratures. Results for the ground-state energy density 
Eq/L for a lattice L = 48 are shown in Fig. [TJ the data 
for Eq/L are essentially converged for this system size 
(within statistical errors, as we can conclude by compar- 
ing against results obtained from smaller lattice sizes, not 
shown here), and good agreement is obtained with the 
results from the LDE and NE described in the previous 
subsection. 



III. NUMERICAL RESULTS 

We present estimates for the location of the quantum 
critical points in the phase diagram of the Hamiltonian 
Eq. ([1]) and results for the excitation spectra and gaps 
obtained by using the numerical procedures described in 
the previous section. 



QMC simulations have been performed by employing 
the ALPS (Algorithms and Libraries for Physics Simu- 
lations) libraries— implementation of the directed loops 
algorith m 40 ! 41 in the Stochastic Series Expansion (SSE) 
framework^ 2 - We have simulated the model Eq. (JTJ) on 
lattices with length L ranging from 24 to 48, applying 
periodic boundary conditions. In order to assess ground- 
state properties, simulations were performed for temper- 
atures T = 1/100 < 1/2L: results for the spin-stiffness ps 
(Fig. EI and the Binder cummulant Q 2 [Eq. IfTT]). Fig. [7] 



A. Large-D — Haldane Phase Transition 

The Gaussian quantum phase transition between the 
large-D and the Haldane phases has been studied in a 
number of recent worksi 20 i 21 ' 24 i 25 The critical point was 
pinpointed by Chen et al^ by performing exact diago- 
nalizations of small clusters with twisted boundary con- 
ditions. Furthermore, these authors verified that the DH 
transition is described by a conformal field-theory with 
central charge c = 1, as further confirmed by Boschi et 
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al£L using a combined ficld-thcorctic and numerical anal- 
ysis. More recently, Tzeng and collaborators; 2 ^ 25 moti- 
vated by the current interest in applying concepts and 
tools from quantum information theory to the study of 
condensed matter systems^ showed that the DH phase 
transition can be located by investigating the scaling be- 
havior of the ground-state fidelity and the entanglement 
entropy, and arrived at the estimate^ D^ H /J = 0.97. 
Here, we investigate the DH transition employing the nu- 
merical methods described in Sec. [H] 

1. Spin Stiffness 

In Fig. Uta) we show our QMC results for the spin- 
stiffness ps, obtained in terms of the winding number w, 
ps = 3(u> 2 )/2/3 (/? is the inverse temperature),— for sys- 
tem sizes L ranging from 24 to 48. Away from the critical 
point, both the large- D and the Haldane phases display 
exponentially decaying spin correlations and therefore we 
expect ps to approach zero in the thermodynamic limit, 
a trend clearly discernible in Fig. [UJa). On the other 
hand, close to the critical point, where the correlation 
length £ diverges, the spin-stiffness is expected^ to scale 
as ps ~ L 2 ~ z ~ d . Since d = 1 is the spatial dimension and 
the dynamic critical exponent is expected to be z = 1 
(see below), one expects ps to assume a size indepen- 
dent value at the transition point (similarly to what hap- 
pens for the 5=1/2 chain with easy-plane anisotropy, 
see Ref. H3) , as confirmed by our QMC data plotted in 
Fig. [UJa). We remark that this "peaked behavior" for ps 
is a signature of a transition between two phases with 
exponentially decaying spin correlations and it should be 
contrasted with the "crossing behavior" for p$L observed 
in more conventional order-disorder quantum transitions 
in two-dimensional systems (see e.g. Refs. l48ll49lj50T ) . 

In order to estimate the location of the quantum crit- 
ical point Z?J=? H /J and the correlation length critical ex- 
ponent v, we assume the scaling ansatz 

Ps (t,L) = f ps (tL 1 ^) , (10) 

with reduced coupling t = (D/Dj? H — 1). Data col- 
lapse is achieved for Dg H /J = -0.971(5) and v = 1.4(1) 
[Fig.J2Tb)], values in good agreement with the ones from 
Ref. [2J (Dg H / J = -°- 97 and v = L42 > 1-45)42. We 
use our result v = 1.4(1) to calculate the Luttinger- 
liquid parameter K by employing the relatio n 21 ' i 25 v = 
1/(2 — K) and obtain K — 1.29(5), consistent with previ- 
ous findings . 21 ' 24 i 25 Additionally, the bulk spin-stiffness 
at the critical point is estimated to be pg H = 1.619(5). 

2. Excitations 

The lowest-lying excitations in the large- D phase lie in 
the iSf = ± 1 sectors and we calculate their dispersion re- 
lation by applying the LDE presented in Sec. Ill A II Re- 
sults obtained after a standard Pade analysis are shown 




9 D/J 



FIG. 3: (Color online) (a) Excitation spectrum in the large- 
D phase, for D/J = —2, obtained from the LDE. (b) Gap 
for excitations (g = n) in the large-D phase as a function 
of D/J, again obtained from the LDE (the series become 
badly behaved beyond the rightmost data point). The con- 
tinuous line has been obtained by fitting the scaling func- 
tion A/J = A [(D° H - D)/J)] z '\ assuming z = 1, to the 
data points depicted as squares, using the QMC estimate 
Dc H /J — —0.971 (vertical dashed line). In both panels, with 
the exception of the rightmost data point in (b), error bars 
are smaller than the depicted symbols. 

in Fig. [3]Ja) for D/J = —2. As we approach the DH 
transition from the large negative D/J side, the gap at 
q = it drops to zero, and the dispersion in that neigh- 
bourhood becomes linear, a trend already noticeable for 
D/J = —2 [Fig. E3a)] and consistent with a dynamic 
critical exponent z = 1. We remark that this behav- 
ior is not reproduced in the strong-coupling analysis by 
Papanicolaou and Spathis^ which gives a quadratic dis- 
persion relation (this is also the case in our results for 
couplings deep into the large-D phase, not shown here). 
The dependence of the gap at q = % on D/ J is shown in 
Fig-Efb). Unfortunately, the series convergence becomes 
irregular beyond the rightmost data point in the figure, 
preventing us from obtaining independent estimates for 
Dq 11 / J and for the critical exponents. In order to par- 
tially circumvent this problem, we fit the scaling func- 
tion A/J = A [(Dg H -D)/J)] zu , assuming z = 1, to the 
data points depicted as squares in Fig. [3^b). By fixing 
D^ H / J = —0.971, as estimated from our QMC data, the 
data is fitted for v = 1.478, a value consistent with the 
QMC result v = 1.4(1) and that further confirms that 
indeed z = 1. 



B. Haldane — Neel Phase Transition 

The Neel phase for large positive values of D/J has 
a twofold degenerate ground-state and therefore the HN 
quantum phase transition is expected to belong to the 
universality class of the two-dimensional Ising model. 
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1. Staggered Magnetization 
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FIG. 4: (Color online) (a) Staggered magnetization per site 
ml as a function of D/J in the Neel phase, obtained from 
the NE by applying a Pade analysis, (b) Location of the 
poles Ac in the Dlog Pade approximants for the staggered 
magnetization series, as a function of D/J. In both panels, 
lines are only guides to the eye. 
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FIG. 5: (Color online) Pseudo-particle representation of the 
S = 1 chain.— Sites with Sf = ±1 are seen as being oc- 
cupied by spin-half particles with pseudo-spin components 
Sf = ±1/2 (depicted as +'s and — 's) and sites with Sf = 
as being empty (occupied by "holes", represented by circles) 
In this language, the Neel ground-state corresponds to an un- 
doped antiferromagnet (a). Excited states are obtained by 
doping the system with holes (b). As depicted in (c), a hole 
doped into the system decays into "spinon" (enclosed by the 
rectangle) and holon (enclosed by the ellipsis) constituents 
[the state depicted in (c) is obtained from the state in (b) by 
allowing the hole to "hop" four sites to the left, as indicated 
by the arrows]. Since a spinon costs an energy ~ J, excited 
states containing only a holon (d) are expected to have lower 
energy (see main text). 



Chen et a/.— located the HN critical point by apply- 
ing a phcnomcnological renormalization group analysis 
to data from exact diagonalization of clusters with up 
to 16 sites. More recently the estimate^ D™ / J = 0.31 
was obtained by Tzeng and Yang^i who studied the scal- 
ing behavior of the ground-state fidelity close to the HN 
critical regime. 



We first employ the NE discussed in Sec. Ill A 21 in 
order to calculate the staggered magnetization, m| = 
{{—l) % Sf)/L, as a function of D/J. The results obtained 
by applying standard Pade approximant extrapolations 
to the series in A are shown in Fig. HJa). The position of 
the HN quantum critical point can be determined by a 
Dlog Pade analysis of as a function of A: as shown in 
Fig.gljb), the estimate D™ / J = 0.295(25) [highlighted 
as the shaded region in Fig. Ufa)] for the HN critical 
point is simply obtained as the range of values for D/J 
consistent with a pole at A = 1, where the full Hamil- 
tonian Eq. (TTJ is recovered. From the Pade extrapola- 
tion we also obtain an estimate^ (3 = 0.147(13) for the 
critical exponent associated to mj, which is somewhat 
larger than, but not incompatible with, the exact result 
(3 = 1/8 for the 2D Ising universality class. We also note 
that there is little sign of ml vanishing in the shaded re- 
gion in Fig. [4|a): this can be explained by the small value 
of the expected critical exponent (3 = 1/8, which implies 
that TOg plunges steeply to zero at the critical point, a 
behavior which naive Pade approximants will hardly pick 
up. 



2. Excitations 

More accurate estimates for Dj? N /J can be obtained 
by analyzing excited states above the Neel ground-state. 
At this point, following den Nijs and Rommelse^! it is 
useful to interpret the S = 1 chain as a diluted system of 
5=1/2 pseudo-particles: as we show in Fig.0 sites with 
Sf = ±1 are seen as being occupied by spin- half particles 
with pseudo-spin components Sf = ±1/2 and sites where 
Sf = as being empty (occupied by "holes" ) . Using this 
language, the Neel ground-state is equivalent to an "un- 
doped" antiferromagnet [Fig. Ufa)] and, for small (pos- 
itive) values of D/J, the low- lying excited states lie in 
the "one-hole sector" (containing one site with Sf = 0). 
Interestingly enough, the situation is reminiscent of spin- 
charge separation in one-dimensional fcrmionic systems 
(see e.g. Ref. [H) where a hole doped into the system 
[as depicted in Fig. \5\b)] fractionalizes into "holon" and 
"spinon" constituents [respectively enclosed by an ellipsis 
and a rectangle in Fig.OJc)]. 

Since one spinon (two consecutive sites with the same 
pseudo-spin component) has an energetic cost ~ J, we 
expect Sf. = states solely displaying a holon to have 
a lower excitation energy, as confirmed by our results 
below. Such a holon state is depicted in Fig. [5fd): we 
replace one pseudo-spin in the Neel configuration shown 
in Fig. [5]Ja) by a hole and flip all pseudo-spins to its 
right. Antiferromagnetic correlations "across the hole" 
allow for the hole to delocalize, lowering the energy of 
holon excited states^ Furthermore, these holon excita- 
tions can be seen as precursors of the ground-state in 
the Haldane phase: as shown by Nijs and Rommclsc^ 
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FIG. 6: (Color online) (a) Spectrum for "holon" excitations 
[Fig. [5jd)] at D/J — 1. The vertical dashed lines indicate 
momenta for which the NE expansion seemingly diverges, (b) 
Spectrum for "holon-spinon" excitations [Fig. [5jb-c)], also for 
D/J — 1. (c) Excitation gaps as a function of D/J. The gap 
for holon excitations (q = it, black circles) vanishes linearly^ 
[u = 1.01(3)] at D% N /J = 0.316(2) as seen in the magnified 
plot (d). Data points in all panels were obtained by applying 
a Pade analysis to the results obtained from the NE. 



VBS states can be interpreted as a fluid with positional 
disorder and long-range antiferromagnetic order. This 
picture is consistent^! with the holon state depicted in 
Fig- 13d), but not with the holon-spinon excitation shown 
in Fig. [5Jb-c)i££ Therefore, the HN critical point can be 
determined by locating the value D™ / J where the gap 
for the holon excitations vanishes. 

Results for the spectrum of holon excitations obtained 
from the NE are shown in Fig. [SJa) for D/J = 1. We 
remark that the scries seemingly diverges at the commen- 
surate momenta q = 2tt/5 and q = 8it/5 [indicated by the 
vertical dashed lines in Fig. [U^a)] and that the results at 
higher energies have a rather poor convergence, as in- 
dicated by the relatively large error bars in the figure. 
This suggests that holon excitations decay into multi- 
holon states and that a continuum of excited states exists 
at high-energies. Unfortunately, the poor convergence of 
the (short) series for two- and three-holon excited states, 
obtained by applying the procedure described in Ref. l54l . 
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FIG. 7: (Color online) (a) Staggered magnetization Binder 
cumulant Q2 [see Eq. for couplings D/J close to the 

Haldane — Neel critical point, as obtained from QMC simu- 
lations for system sizes L ranging from 24 to 48. (b) Data 
collapse is achieved for D% N /J = 0.312(6) and v = 1.1(1). 



prevents us from further analyzing this issue, which is 
left open to future investigations. On the other hand, 
our results for low-energy excitations around q = tt are 
nicely converged, allowing us to precisely locate the crit- 
ical point. In Figs.[BJc-d) we show the dependence of the 
gap A/ J for holon excitations on D/J and we can see it 
vanishes at D™ / J = 0.316(2), with a critical exponent— 
v = 1.01(3) (directly obtained from the Pade analysis) 
consistent with the 2D Ising universality class. 

For the sake of comparison, we also show the disper- 
sion relation for the Sj, = ±1 holon-spinon [Fig. E^b- 
c)] in Fig. [B^b), also for D/J = 1. The dependence of 
the gap at q = on D/ J is shown in Fig. [6tc) . Note 
that our results indicate that the energy of this excita- 
tion remains finite at the transition, confirming that a 
holon-spinon state has a higher excitation energy than 
an isolated holon state. As we see in the plot Fig. [6](c) , 
for large values of D/J the difference in energy between 
holon-spinon and holon excited states is approximately 
J, confirming our qualitative analysis above. 



3. Binder Cumulant 

The previous findings arc confirmed by our QMC re- 
sults for the second order Binder cumulant Q2 for the 
staggered magnetization mj, given by 



(K) 



z\4\ 



3<(mf) 



2\2 



(11) 



Q2 is expected to display universal behavior in the critical 
regime and results obtained from different lattice sizes 
should cross close to the critical point, as confirmed by 
the data shown in Fig. [7{ a )- In order to estimate the 
critical point we assume the scaling ansatz Q2(t, L) = 
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fq^tL 1 ^), with t = {D/D™ - 1). Data collapse is 
achieved for D™ / J = 0.312(6) and v = 1.1(1), as shown 
in Fig. El^b), consistent with the estimates obtained from 
the NE. A slightly smaller value, D^/J = 0.310(6), is 
estimated from a data collapse obtained by fixing the 
critical exponent to the known value for the 2D Ising 
universality class, v = 1. 

Additionally, from the data collapse displayed in 
Fig. [7](b) we arrive at the estimate Q* 2 = 0.515(10) for 
the value assumed by the Binder cummulant at the crit- 
ical point, markedly different from the known result^ 
Q2 = 0.61069... for the two-dimensional Ising universal- 
ity class. This discrepancy can be ascribed to differences 
in the way Q2 is calculated in classical and quantum 
Monte Carlo simulations: although the d-dimensional 
quantum system is indeed formally equivalent to a clas- 
sical model in d + 1 dimensions, equal-time expecta- 
tion values for diagonal operators are averaged along 
the imaginary-time direction in quantum simulations. In 
other words, our results for should be consistent with 
the ones obtained for the classical two-dimensional Ising 
model if the moments of the magnetization appearing 
in the definition Eq. ([lip are evaluated along individ- 
ual lines on the square lattice. Since, to the best of our 
knowledge, no such calculation has been done, this pre- 
vents us from comparing our estimate Q\ = 0.515(10) 
with published results. We remark that a similar situ- 
ation was observed in Ref. \4g for S = 1/2 bilayers, ex- 
pected to belong to the universality class of the classical 
Heisenberg model in three dimensions. 



C. Haldane Gap 

Finally, we have investigated the SFE introduced by 
Singh^ and discussed in Sec. Ill A 31 in order to estimate 
the Haldane gap at D/J = 0. Singh assumes that the 
magnon gap in the Neel phase of the SFE system will 
extrapolate continuously to the Haldane gap as A — ► 1, 
in spite of the phase transition that must occur at or 
before A = 1. The extrapolated value for the magnon 
gap, A/ J = 0.42(1)^. obtained from a series involving 
terms of up to A 10 as presented in Ref. HH, is indeed 
consistent with the high-precision result for the Haldane 
gap A/ J = 0.41050(2) from density-matrix renormal- 
ization group (DMRG) calculations^ However, higher- 
order Pade approximants for an extended series compris- 
ing terms up to A show a clear trend towards smaller 
values and we arrive at the considerably lower value 
A/ J = 0.23(6). It is thus obvious that a more careful 
analysis is required. If we apply a change of variable 
<5 = 1 — \/l — A to the series, following Huse,— then Pade 
approximants for the series in 5 show a vanishing magnon 
gap, A/ J = 0.00(1), at the transition point A = 1. The 
failure to yield correct results for the Haldane gap, com- 
bined with the fact that it systematically underestimates 
the ground-state energy towards the DH critical point, 
as mentioned in Sec. Ill A 31 in connection with Fig. [TJ 



strongly suggests that the SFE approach is not appro- 
priate for estimating results for the integer-spin Haldane 
systems. 

We recall that a Neel state is chosen as the unper- 
turbed state when performing the SFE, and the VBS 
character of the Haldane phase is ignored. We expect 
physical results to be obtainable from a different expan- 
sion assuming the AKLT model [/3 = 1/3 in Eq. ^] as 
the unperturbed Hamiltonian and starting from the exact 
VBS ground-state^ at /? = 1/3, treating A = (/3 - 1/3) 
as the expansion parameter. Such an expansion would 
require the application of the linked-cluster formalism to 
valence-bond states, something that, to the best of our 
knowledge, has not yet been tried and that would consti- 
tute an interesting topic for future work. 



IV. CONCLUSIONS 

Summarizing, we have investigated the one- 
dimensional S = 1 antiferromagnet with single-ion 
anisotropy, described by the Hamiltonian Eq. |TJ, by 
means of linked-cluster series expansions and QMC 
simulations. Our estimates for the zero-temperature 
phase transitions in this model are more precise than 
previous ones and could be used as benchmarks in future 
explorations of the applicability of quantum information 
tools to the study of quantum critical phenomena. 

Our best estimate for the DH critical point [-D§ H / J = 
0.971(5)] has been obtained from a scaling analysis of 
the spin-stiffness ps from QMC simulations. The spin- 
stiffness remains finite at the transition in the thermody- 
namic limit [ps H = 1.619(5)] and vanishes elsewhere, im- 
plying a "peaked behavior" for finite systems. Our result 
for £>q H agrees with the estimate^ obtained by Tzeng et 
al.— [-D§ H = 0.97], the same being true for the estimate 
for the correlation length critical exponent: v = 1.4(1) 
(present work) and v = 1.42, 1.45 (Ref. [H). We have 
further obtained results for the dispersion relation and 
gap for excitations in the largc-Z? region that may be 
of direct experimental relevance given, for instance, the 
current interest on the large- D compound DTN . 11 ' 12 Our 
series for the excited states in the large- D phase have a 
considerably more extended range of applicability com- 
pared to the strong coupling results of Papanicolaou and 
Spathis^ and are available on request. 

Precise results for the HN phase transition have been 
obtained from a linked-cluster expansion (NE) for holon- 
like excitations. By analyzing the gap for these ex- 
cited states we have arrived at the estimates D™ / J = 
0.316(2) and v = 1.01(3). The former compares well with 
the result Dg N / J = 0.31 from Ref. H3 while the latter 
confirms that the HN transition belongs to the univer- 
sality class of the two-dimensional Ising model, with the 
exact exponent v = 1. 

Finally, we have shown evidence that the SFE expan- 
sion of Singh22 does not converge to the Haldane gap at 
D/J = 0. Instead, we propose an expansion around the 



8 



AKLT state [the exact ground-state for the Hamiltonian work has been supported by the Australian Research 
Eq. @ with (3 = 1/3] j2S which we hope to explore in Council, 
future work. 
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